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1 Introduction 



Discrepancies are measures of non- uniformity of point sets that play an important role 
in the Quasi- Monte Carlo method of numerical integration [|I|. A certain class of these 
discrepancies, the so called quadratic discrepancies, can be defined as an average-case 
complexity over a class of functions, the problem class ^, ^. In a number of publications 

1^, 1^ the problem of calculating the probability distribution of discrepancies for sets of 
N truly random points has been tackled. One of the results was the calculation of the 
asymptotic distributions in the limit of infinite N. 

In [H, we introduced techniques from quantum field theory (QFT) to calculate the 
moment generating function G of the probability distribution, suitable to calculate it as 
a series expansion in In this paper, we extend the formalism by the introduction 

of fermions, and give the explicit diagrammatic expansion of logG up to and including 
0{1/N'^). For the Lego discrepancy, the L2-discrepancy in one dimension and the Fourier 
diaphony in one dimension, we give the explicit correction. 



2 The formalism 

We start with a short repetition of the formalism of and continue with the introduction 
of some new tools. 



2.1 Quadratic discrepancies and path integrals 

We shall always take the integration region to be the s-dimensional unit hypercube K = 
[0, 1)"*. The point set X^r consists of points e K, = 1, 2, . . . , A^. Defined as an 
average-case complexity on the problem class of functions : K i— > R with measure /x, a 
discrepancy Z^tv of the point set X^r is given by 

f 1 ^ f 

Dm = N j r/^[0]d/i[0] , r/^[0] = j^^<Pi.^k) - j <P{x) dx . (1) 

k=i 

It is the quadratic integration error, averaged over the problem class. The probability 
density H of the discrepancy as a random variable of random point sets is calculated as 
the inverse Laplace transform of the moment generating function G: 



+ 100 

H(DN = t) = — f e-''G{z)dz , G{z) = ETe"^^] 
2ni J 



(2) 



Here, E denotes the expectation value of a random variable. 

The measure fi is assumed to be Gaussian and in Q it is shown that the generating 
function is therefore given by 

G{z) = j (^j e^['^(")-/K*(f)'^?/]c^a;^ dfi[(P] , g = . (3) 
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In 1^ we suggested to put the A^**^ power in the exponential and interpret G{z) as an 
Euchdean path integral 

G{z) = I V<j>exp{-S[<P]) (4) 



with an action 



S[(f)] = ^ [ (j){x)A{x,y)(f){y)dxdy- N log ( [ e^^'^'^^'^-^^'^'^y^'^yUx] , (5) 



K2 \JK 



where A is the symmetric linear operator with boundary conditions which is the inverse of 
the two-point Green function under the measure fi: 

^{.Xi,y)C{y,X2)dy = 6{xi - X2) , C{y,X2) = / 0(?/)0(a;2) . (6) 



K 



The "infinitesimal volume element" V(f) can be seen as defined by the rule that dfi[(j)] = 
Vcf) exp(— S'o), where 5*0 is the action with z = 0. For notational convenience we put T>(j) 
to the left of the exponential. 

One of the features of this formalism is that the action has a gauge freedom] a global 
translation Oc : (p{x) 1— *• (f){x) + c, c G R leaves ?7Ar[0] and a fortiori G{z) invariant, and 
results in a change of the action at most linear in 0: 

S[0,4>] = S[4>] + acx[4>] + with + H = x[<Pi] + x[H • (7) 

As a result of this, the path integral can be generally written as 

G{z) = y^y^ j V<P eM-nm) " Se[<P]) , Se[<P] = S[<t>] - , (8) 

where F is restricted such that 

POO 

I[F] = / exp(-F(c))rfc (9) 



exists, and ^ is only restricted such that ^[Oc(p] = C,[4>] + c. It is, for example, possible to 
take = j^(f){x) dx and F{x) = Mx^ with M —* 00. In this gauge, will vanish 
from the action, which is reflected in a two-point function that integrates to zero with 
respect to each of its variables. We shall refer to this gauge as the Landau gauge. 

Clearly, the first equation of (^) cannot be satisfied in the Landau gauge, because then 
C integrates to zero Q. To see what happens, we assume that the problem class is a vector 
space with a basis of at least square integrable functions so that a member of the 

problem class can be written as 

= y^0nMn(x) , 0n € R , (10) 



^There is a misprint in Eq. (73) of [Q with respect to this, where S{x — y) should be replaced by 
<5(a; - y) - 1- 
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and that the Gaussian measure on this function space is defined by 



rf/^M = n '"P^ ^^-^ d0„ , a. e R . (11) 

For the measure to be suitably defined, the strengths (T„ have to satisfy certain restrictions 
which can be translated into the requirement that EfD^r] exists. They are the inverse of 
the eigenvalues of C and the basis consists of the eigenfunctions. Therefore, C and A can 
be expressed in terms of the basis: 



where the last equation holds if the basis is orthonormal. The boundary conditions satisfied 
by C and A are those satisfied by the basis functions. With different gauges come different 
bases and strengths. We call a gauge in which the basis is orthonormal a Feynman gauge. 
If the Landau gauge is used, in which Jj^ 0(a;) dx = 0, then the basis functions have to 
integrate to zero: 

[ u^^\x)dx = Vn , (13) 

where the label L indicates the Landau gauge. This means that the basis cannot be 
"complete" in the sense that '^n'^n{^)'^n{y) = ^{^ — u), but that we have 

J]4^)(a;)4^)(|/) = S{x-y)-l . (14) 

n 

The zero mode is isolated and integrated out of the path integral. We want to stress that 
the gauge freedom is something that comes from the original measure /i, and that the 

Landau gauge exists for every quadratic discrepancy. It is a result of the fact that an 
integration error is the same for integrands that differ only by a constant. 

From now on we will denote the two-point function in the Landau gauge by B. It 
satisfies 

f A^''\x,,y)B{y,X2)dy = 5{x,-X2)-l , (15) 

and the discrepancy can then be written as 

1 ^ 

Dn = . (16) 



fcj=i 
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2.2 Fermions as tools to calculate the l/AT corrections 

In m we suggested to make a straight forward expansion in of exp(— 5) to calculate 
G perturbatively. This way, however, the calculation of the perturbation series becomes 
very cumbersome, and the reason for this is the following. We want to use the fact that an 
expansion in corresponds to an expansion around = of the part of the action that 
is non-quadratic in (p. The subsequent terms in the expansions are therefore proportional 
to moments of a Gaussian measure, and can be calculated using diagrams (cf. ||10|)- These 
diagrams, the Feynman diagrams^ consist of lines representing two-point functions and 
vertices representing convolutions of two-point functions. Because the action is non-local, 
i.e. it cannot be written as a single integral over a Lagrangian density because of the 
logarithm in Eq. (j^), the total path integral, thus the total sum of all diagrams, cannot be 
seen as the exponential of all connected diagrams, and it is this that makes the calculations 
difficult. 

In order to circumvent this obstacle, we introduce 2N Grassmann variables T/ij and 
'(/'j, z = 1, 2, . . . , A^. They all anti-commute with each other and commute with complex 
numbers: 

V^iV'j V'jV'i = , ipiipj + ipjifji = , ^pi^pj + 2pj2pi = i,j = l,2,...,N (17) 
cipi - ipiC = , cipi - ipiC = i = 1,2, . . . ,N , c G C . (18) 

Now we use the well known Gaussian integration rules for Grassmann variables to write 
the N^^ power in Eq. (^ as an exponent and get 

G{z) = j V(PDipDiP exp{-S[(l),'ip,ip]) , (19) 

TV 

S[^,i>,ij] = 11 cP{x)A{x,y)(P{y)dxdy+ I e^l'^^^^-^K-^fe)'^?/] rf^; V]^^^^ , (20) 
2 Jk2 Jk ^ 

where we introduced the notational shorthand 

DTp = dTpidip2 ■ ■ ■ dipN , Dip = dipidip2 ■ ■ ■ dipN ■ (21) 

Notice that this action contains the same gauge freedom, so that the action becomes 
completely local if the Landau gauge is used. We have to introduce the "fermion fields" to 
achieve this, but for calculating the perturbation expansion they are much easier to handle 
than the logarithmic potential. From the action (^) we obtain the Feynman rules. To 
calculate a term in the 1/A^-expansion of G, the contribution of all diagrams that can be 
drawn using the Feynman rules and carry the right power of 1/A^ has to be calculated. 

Because in this paper we want to calculate the path integral itself, and no correlation 
functions, we only have to consider vacuum diagrams, i.e., diagrams without external 
legs. Furthermore, we will always use the Landau gauge, because then action is local, so 
that we only need to calculate connected diagrams. The whole generating function is the 
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exponential of the sum of these diagrams. The contribution of the connected diagrams we 
denote by W, so that 



G{z) 



(22) 



The diagrammatic expansion is an expansion in the subsequent terms of which we 
denote by 



W{z) = Wo{z) + ^W,iz) + ^W,iz) + ■ ■ 
In the Landau gauge the rules can be summarized by 



boson propagator: x ■ 
fermion propagator: i 



y = B{x,y) ; 
■j = 5, 



vertices: 



-g^ X convolution , p > 2 



(23) 



(24) 
(25) 

(26) 



In the vertices, boson propagators are convoluted as f-^B{y,Xi)B{y,X2) ■ ■ ■ B{y,Xp) dy, 
fermion propagators as J2f=i^ii,j^j,i2^ then these convolutions are multiplied. As a 
result of this, the bosonic part of each diagram decouples completely from the fermionic 
part. The contribution of the fermionic part can easily be determined, for every fermion 
loop only gives a factor — A^. The main problem is now to calculate the remaining bosonic 
part. Finally notice that, as a result of the Landau gauge, vertices with only one bosonic 
leg do not exist. 



3 The diagrammatic expansion 

We want to stress again that we only need to calculate the connected diagrams. The sum 
of the contributions of all these diagrams gives W = logG. Usually, a Feynman diagram 
is a mnemonic representing a certain contribution to a term in a series expansion, i.e. a 
label. We will use the same drawing for the contribution itself, apart of the symmetry 
factor of the diagram. For example, the contribution of the diagram ( )(y is equal to 
^Ng'^ B{x, x) dx and its symmetry factor is equal to |, so that we write 

^C/;0 = ^ iB{x,x)dx . (27) 



K 



^In y we called another two-point function the propagator, namely the inverse of A(x, y) — 2z6{x^y)+2z, 
denoted by Qz{x,y). This one is obtained from B{x,y) by the re-summation of certain diagrams, as we 
will show in the following. 
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3.1 The zeroth order 



The contribution to the zeroth order in can only come from diagrams in which the 
power of coming from the vertices cancels the power of A^, coming from the fermion 
loops. This only happens in diagrams with vertices with two bosonic legs only, and in which 
the fermion lines begin and end on the same vertex. To write down their contribution, we 
introduce the two-point functions Bp, p = 1,2, ... , defined by 

Bi{xi,X2) = B{xi,X2) , Bp+i{xi,X2) = / Bp{xi,y,)B{y,X2)dy . (28) 

Jk 

The zeroth order term is given by 





i(";;o + iO;;;;:o + lA'jt^ + ■■■ = f B,(.,.)d. . (29) 



The factor l/2p is the symmetry factor of this type of diagram with p fermion "leaves". 
If we substitute g = ^j2zjN in this expression, we find exactly the result of Eq. (21) in 
§]. If we use the spectral representation of B and assume that the basis functions are 
orthonormal in the Landau gauge, we get 

^o(^) = -i5^1og(l-2za^J , (30) 

n 

where L indicates the Landau gauge. This expression is the same as Eq. (68) in 0. 
3.2 The first order 

As we have seen before, bosonic two-point vertices with a closed single fermion line con- 
tribute with a factor 2z, and without any dependence on A^. Therefore, it is useful to 
introduce the following effective vertex 



,....0. 



1 p = 1 ^ p = ^9 X convolution , (31) 

-''''1 / .''''I 
2' / ; 2' / ; 

4 4 

and the following dressed boson propagator 

X y = X y + X ♦ y + x y 

+ X y + ■ ■ ■ 

00 

= Y.{2zr'Bp{x,y) . (32) 
p=i 



6 



If we assume that the basis is orthonormal in the Landau gauge, we can write 

2 

X y = E l ^9'% ^^n\^Wn\y) , (33) 

which is, apart of a factor 2z, the same expression as in Eq. (67) in |p. The dressed 
propagator is equal to the propagator in the Landau gauge as we defined it in Eq. (24) in 
i: 

X y = gi''\x,y) . (34) 

This is easy to see, because it satisfies 



xi - X2) - I , (35) 



„ 00 

/ [A^'^\x,,y)-2z6{x,-y) + 2z]y2i^zy''^pi^^y)dy = S{ 

p=i 

just hke ^i^'' by definition. Notice that ^l'"^ and B satisfy the relation 

\iuig^^\x,y) = Gf2o{x,y) = B{x,y) Wx,y,eK . (36) 

Furthermore, notice that Qi^'^ and Wo satisfy 



^ Woiz) = I gf\x,x)dx , (37) 

K 



dz 



and that this relation determines Wq uniquely, because we know that M/o(0) has to be equal 
to in order for the asymptotic probability distribution to be normalized to 1. From now 
on, we will omit the label L. 

The first order term in the expansion of W{z) is 

>'.W-ioO + ioOO + i@nCX5 + ^0. (38) 

or, more explicitly, 

Wi{z) = yJ '3z{x,xfdx - Y (^j Gzix,x)dx^ - z^ j ^g^{x,yf dxdy 

+ z^ [ g^{x,x)Q^{x,y)g^{y,y)dxdy + ^ /" g^ix^yfdxdy . (39) 

3.3 The second order 

The second order term in the expansion of W{z) is denoted by ■^W2{z) and is given by 
1 U 1 




48 CyVj 48 





+ 
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+ 



+ 




+ sCxX>o + i!5 






1 

+ 4 





+ 



+ 






3.4 One-vertex decomposability 

For some discrepancies, the contribution of a bosonic part of a diagram that consists of two 
pieces connected by only one vertex, is equal to the product of the contribution of those 
pieces. Such diagrams we call one-vertex reducible, and discrepancies with this property 
we call one-vertex decomposable. Examples of such discrepancies are those for which E 
is translation invariant, i.e., B{x,y) = B{x + a,y + a) \/x,y,a G K, such as the Fourier 
diaphony. Also the Lego discrepancy with equal bins is one-vertex decomposable. In 
contrast, the Lg-discrepancy is not one- vertex decomposable. 

As a result of the one-vertex decomposability, many diagrams cancel or give zero. For 
example, the first and the second diagram in cancel, and the fourth gives zero, so that 




(40) 
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To second order, only the following remains: 



Ar2 



+ 



+ 





(41) 



We now derive a general rule of diagram cancellation. First, we extend the notion of 
one- vertex reducibility to complete diagrams, including the fermionic part, with the rule 
that the two pieces both must contain a bosonic part. Consider the following diagram 



0<K) 



(42) 



The only restriction we put one the "leave" A is that it must be one- vertex irreducible with 
respect to the vertex that connects it to the fermion loop. For the rest, it may be anything. 
We define the contribution of the leave by the contribution of the whole diagram divided 
by — A^, and denote it with C{A). This contribution includes internal symmetry factors. 
Now consider a diagram consisting of a fermion loop as in diagram (|4^) with attached to 
the one vertex rii leaves of type Ai, n2 leaves of type A2, and so on, up to Up leaves of type 
Ap. The extra symmetry factor of such a diagram is {ni\n2\ ■ ■ -np!)"^, and, for one-vertex 
decomposable discrepancies, the contribution is equal to the product of the contributions 
of the leaves, so that the total contribution is given by 



q=l 



C{A, 



n 



(43) 



q- 



Now we sum the contribution of all possible diagrams of this kind that can made with the 
p leaves, and denote the result by 



O 



^ E n 



C{A, 



^1 1 ^q- 

ni,n2,...>l <?=1 



(44) 



Because the black square in l.h.s. of Eq. (^^ represents all possible ways to put the leaves 
together onto one vertex, the sum of all possible ways to put the leaves onto one fermion 
loop is given by 



n=l 



n 



(45) 



The (—1)" ^ in the sum comes from the vertices and 1/n is the extra symmetry factor of 
such diagram with n vertices. The sum can be evaluated further and is equal to 

p 

-iVloge^?=i^(^') = -NY^C{Aq) , (46) 

9=1 

i.e., the sum of all possible ways to put p different leaves onto one fermion loop is equal 
to the sum of all leaves, each of them put onto its own fermion loop. This means that 
diagrams, consisting of two or more leaves put onto one fermion loop, cancel. 

Now consider the following equation, which holds for every one-vertex decomposable 
discrepancy: 




where we only assume that B is not of the type on the l.h.s. of Eq. (^5]). The minus 
sign comes from the fact that the first diagram has one vertex less. Because the number 
of fermion lines, a fermion loop consists of is equal to the number of vertices it contains, 
we can always pair the diagrams into one diagram of the l.h.s. type and one of the r.h.s. 
type so that they cancel. We can summarize the result with the rule that for one-vertex 
decomposable discrepancies, only the one-vertex irreducible diagrams contribute. 



4 Applications 

We apply the general formulae given above to the Lego discrepancy, the Lg-discrepancy in 
one dimension and the Fourier diaphony in one dimension. 



4.1 The Lego discrepancy 

For the Lego problem class, for example defined in the basis consists of a set of char- 
acteristic functions ^n, n = 1,2, . . . ,M of M disjunct subspaces of K. We will denote the 
measure Jj^ '^n{x) dx of subspace n by w„ and we have X^nli = 1- The strengths (j„ are 
equal to 1/ ^Jw^ , and the propagator is given by 

^ ^^^M£RM_i. (48) 

n=l 

With this choice of (T„, the discrepancy is just the statistic that determines how well the 
points are distributed over the bins. It is easy to see that Bp{x, y) = B{x, y), p = 2,3, . . . , 
so that the dressed propagator is given by 

Qz{x,y) = -^13{x,y) . (49) 



10 



In [Q, the propagator is given as an M x M- matrix with matrix elements 

5r, 



l-2z 



Wr. 



- 1 



(50) 



and acting in the M-dimensional function space, rather than as a two-point function. This 
follows naturally from the path integral, which is an M-dimensional integral. The obvious 
and correct relation between the two is that 



M 



g.{x,y) = Ux)gi%^m{y) ■ (51) 

n,m=l 

The zeroth order term can be found with the relation of Eq. (p7D, which results in the 
following expression 

1 f M - 1 

Wo{z) = --log(l-2z) J^Bix,x)dx = -_log(l-2z) , (52) 

in agreement with Eq. (44) in To write down the first order term, we introduce 

^ 1 2z 
M2 = ^ — , and vi^) = —— , (53) 

so that 



Wn L — 2z ' 

n=l " 



Wi{z) = - (Ma - - 2M + 2) 7]{zf + — (5M2 - 3M^ - 6M + 4) r]{zf . (54) 
8 24 

If the bins are equal, so that Wn = 1/M n = 1, 2, . . . , M, then only the contribution of the 
diagrams of Eq. (^^ remains, and the result is 

W^iz) = -^Evizf + ^iE'-EUzf , (55) 

where we denote 

E = M -1 . (56) 
To second order in 1/A^, the contribution comes from the diagrams in Eq. (^l]) , and is given 

by 

W2{z) ={5E^-12E^ + 7E)^^ + {E' - 6E^ + BE)"^ 

48 8 

+ (e3_28e2^43)^ + {-E'-BE)"^ . (57) 

In Appendix A, we present the expansion of G{z) in the case of equal bins, up to and 
including the 1/A^^ term. It is calculated using the path integral expression (^) of G{z) 
and computer algebra. The reader may check that this expression for G{z) and the above 
terms of W{z) satisfy G{z) = e'^*-^-' up to the order of l/N"^. 
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4.2 The L2-discrepancy 

For the Lj-discrepancy in one dimension, for example defined in the gauge freedom 
is a freedom in the boundary conditions which the members of the problem class have to 
satisfy. A acts on the members as 

(/I0)(x) = -^(^) , (58) 
and in the Landau gauge, the boundary conditions are given by 

{x)dx = , ^(0) = ^(1) =0 . (59) 



L 



The basis in the Landau gauge is given by the set of eigenfunctions of A with the boundary 
conditions above, which is {v^ cos(r;,7rx), n = 1, 2, . . . }, so that the propagator is given by 



^, , v-^ ^ l-Obl 74/1 X I CObI 74/1 O I , , 19 19 1 / s 

B{x,y) = = mm(x,2/)-x + |a;2-j/ + iy2+ 1 _ (gg) 

n=l 

The dressed propagator is given by 



^ 2 cos{mTx) cos{mTy) 
yz[x,y) = > ^ 61 

n=l 

= -^~7r—- {cos[u{l - \x + y\)] + cos[u{l - \x - y\)]} , (62) 

u sin u 

with 

u = . (63) 
The zeroth order term can be obtained using Eq. (|37|): 



Woiz) = -llogf^^) , (64) 



2 \ u 

which is the well-known result. After some algebra, also the first order term follows: 

W,{z) = ^^24-8^-7^-7^-2^') . (65) 
288 \ smu sin M tanw tair u J 

4.3 The Fourier diaphony 

Usually, the Fourier diaphony is defined in terms of a basis that is in the Landau gauge 
already. It is, for example, given in |^, and in one dimension, the basis is given by the 
functions 

'U2„_i(x) = v^sin(27rnx) , M2n(a;) = V^cos(27rnx) , n = l,2,... , (66) 
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and for the strengths we take 



a2n-l = (T2n = - , U = 1,2, . . . . (67) 

n 

Notice that the basis functions satisfy —-^Un{x) = ^Un{x), so that, from this point of 
view, the only relevant difference between the L2-discrepancy and the Fourier diaphony 
are the boundary conditions on the members of the problem class. 
The propagator is given by 

,) = ± 2cos(2,,.{x-y}) ^ _ gf, _ _ _ ^j,, _ (gg, 

Tl O 

n=l 



where we use the notation {x} = x mod 1. The dressed propagator is given by 

r ,,\ - V" 2cos(2r27r{a:-y}) _ tt^ / v cos[v{2{x - y} - 

yA^.y) - 2^ n2-2^ ~ vA 2sint; i ' ^ ' 

n=l ^ ' 

where 

V = V27i^z . (70) 

This two-point function is, apart of a factor vr^/f^, the same as the one in Eq. (26) in [Q. 
The zeroth order term can easily be obtained from the dressed propagator and is given by 

Woiz) = -log(^) , (71) 

which is in correspondence with Eq. (21) in 0. Because the propagator is translation 
invariant, i.e., B{x + a,y + a) = B{x,y) Wx,y,a G K, the contributions of the first two 
diagrams in Eq. ( p8|) cancel, and the contribution of the fourth diagram is zero. The 
contribution of the remaining diagrams gives 

W,{z) = ^(s + v'-S^) . (72) 
6b \ sm V I 



5 Conclusions 

In addition to the machinery of QFT introduced in , we introduced fermions to calculate 
the moment generating function G of the probability distribution under sets of random 
points of a quadratic discrepancy Da?. They allow for an expansion in the inverse of the 
number of points of the logarithm W of G, where the contribution to each term in the 
expansion can be represented by a finite number of connected Feynman diagrams. We 
have presented the diagrams up to the order of 1/A^^ for the general case, and derived a 
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rule of diagram cancellation in the case of special discrepancies, which we call one-vertex 
decomposable. 

We have applied the formalism to the Lego discrepancy, the Lg-discrepancy in one 
dimension and the Fourier diaphony in one dimension, and calculated the first two terms 
Wo and Wi/N in the expansion. For the Lego discrepancy, this resulted in Eq. (|5^ ) and 
Eq. (0), for the Lg-discrepancy in Eq. (^) and Eq. (|65|), and for the Fourier diaphony in 
Eq. (^) and Eq. (|72D . The Fourier diaphony and the Lego discrepancy with equal binning 
are one- vertex decomposable. For the latter, we also calculated the term W2/N'^, which 
is in correspondence with the result of an alternative calculation up to the order of 
given in Appendix A. 

Calculations become very cumbersome for high orders because of the number of dia- 
grams involved. A situation in which the formalism can still be powerful is when another 
parameter in the definition of the discrepancy, such as the dimension of the integration 
region or the number of bins in case of the Lego discrepancy, becomes large. This param- 
eter can then serve as an extra order parameter in the determination of the importance of 
the contribution of the diagrams, and can lead to a substantial reduction in the number 
of relevant diagrams. In 0, we will present our results with respect to this for the Lego 
discrepancy. 



Appendix A 

If we define, for the Lego-discrepancy with equal bins, E = M — 1, rj^z) = 2z/{l — 2z) and 

(1 - 2.)^/^G(.) = J2^^n\E) , (73) 

?i,p>0 

then the only non-zero Cn\E) up to = 4 are given by 

Cf\E)=E(-E'--E + -\ 
^ ^ ' V48 48 Sy 

Cf\E) =e( —E' + -E' -J^E + -] 
^ ^ ' V288 72 288 48/ 

Ct\E)=E(~^E^ + ^E~^-^\ 
^ ^ ' V 48 12 48y 
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3 ^ ^ V240 30 120 240 / 

ci'\E) = Eh^E'-'-l^E^+'-^E-'-^) 

^ ^ ' \576 384 576 48 / 

^ ^ ^ V576 1152 1152 576 48 J 

^(8W^x ^/^ 11 ^4 85 o 5125 . 1555 ^ 17\ 

" ^ ^10368 3456 3456 10368 576 72 

_ 1 , 37^2 337„ 1397 



cf\E) = E^E- - 515b^. + Z155b^ _ 15^55^ + 51551) 

^ ^ ^ Vl440 960 720 320 1440 J 

cP(E) = e(^e' - H5555£. + 5Z!152b^ _ 511151s + 51555) 

^ ^ ^ V960 5760 5760 2880 144 / 

cf (£) = e(^e^ + 155Z5£. _ 5555551^3 + 1446743^, 

^ ^ ^ V 23040 23040 92160 7680 

9583187^ 294773 \ 

ip _|_ I 

23040 1152 J 

^(9),^, 73 . 35077 . 781079 . 993515 . 564301 ^ 24607 

Cr(E) = E{ + E^ E^ + E^ E + 

^ ^ ^ V6912 13824 13824 3456 1152 96 

cr(E) = e(^e^ + ^E^ + 155151E* - 55515Zb3 

^ ^ ^ V 13824 3072 34560 9216 

1653251 . 253799^ 145199\ 
H E^ E + 



6912 768 960 



cl")(B) = e( J^e' + ^E^ + 55555s4 _ 1555515^3 

4 ^ ^ V 41472 13824 13824 41472 

19783 . 137875 ^ 1565 

H E^ E H 

192 1152 32 

^/ i p. 11 2431 ^5 . 155735 



Cr'(E) = E{ E' + E^ + E^ + E"" 

^ ^ ' V497664 31104 82944 124416 

3942431 . 249239 . 250141 ^ 2575\ 

497664 13824 13824 384 J 
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